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ABSTRACT 

We examine the expected signal from annihilation events in realistic cold dark mat- 
ter halos. If the WIMP is a neutralino, with an annihilation cross-section predicted 
in minimal SUSY models for the lightest stable relic particle, the central cusps and 
dense substructure seen in simulated halos may produce a substantial flux of energetic 
gamma rays. We derive expressions for the relative flux from such events in simple 
halos with various density profiles, and use these to calculate the relative flux pro- 
duced within a large volume as a function of redshift. This flux peaks when the first 
halos collapse, but then declines as small halos merge into larger systems of lower 
density. Simulations show that halos contain a substantial amount of dense substruc- 
ture, left over from the incomplete disruption of smaller halos as they merge together. 
We calculate the contribution to the flux due to this substructure, and show that it 
can increase the annihilation signal substantially. Overall, the present-day flux from 
annihilation events may be an order of magnitude larger than predicted by previ- 
ous calculations. We discuss the implications of these results for current and future 
gamma-ray experiments. 

Key words: elementary particles - dark matter - galaxies: structure - gamma rays: 
observations - gamma rays: theory 



1 INTRODUCTION 

Dark matter is omnipresent in the universe. Most of it is 
non-baryonic, and a favoured candidate is a weakly inter- 
acting massive particle (WIMP), often generically taken to 
be the lightest stable relic particle surviving from when the 
universe was supersymmetric. The freeze-out of such a neu- 
tralino of mass m x occurs at kT ~ m x /20, and the annihila- 
tion cross-section determines the current value of the CDM 
density Q, c dmh 2 . For minimal SUSY, for example, one can 
derive the relation between annihilation cross-section and 
particle mass, and compute the various annihilation prod- 
ucts, including continuum gamma rays from tt° decays and 
line gamma rays from rare quark decays (c.f. Bergstrom 2000 
for a recent review). While cosmological observations spec- 



ify Q c 



0.1, scanning over minimal SUSY parameter 



space results in an uncertainty in the gamma-ray emissivity 
of several orders of magnitude. 

Our galactic halo is a logical place to look for evidence of 
annihilations. Unfortunately, for a uniform dark halo with a 
realistic density profile, even the most optimistic models fall 
short of the observed diffuse high-galactic-latitude gamma- 
ray flux, as measured by EGRET, by an order of magnitude 
or more (Ullio et al. 2002). In fact, high-resolution numer- 
ical simulations show that the dark halo has considerable 



substructure (Klypin et al. 1999; Moore et al. 1999). This 
substructure may boost the annihilation flux substantially. If 
this is indeed the case, then the isotropic diffuse background 
flux from the many small halos that merged in the past into 
our halo and others may also become significant. There is 
considerable uncertainty among cosmological halo simula- 
tors, however, about the quantitative role of substructure 
in our dark halo, and of the concentration of the substruc- 
ture and of the dark matter itself towards the centre of the 
galaxy. 

No reliable estimates have been given up till now of the 
properties of halo substructure on very small scales, and in 
particular of their dependence on halo mass and their evo- 
lution with cosmological epoch. We have developed a semi- 
analytical model of halo formation which is capable of fol- 
lowing the key physics of tidal disruption of substructure 
during merging. In principle this approach has arbitrarily 
high resolution. Hence we are able to provide robust cal- 
culations of the annihilation flux generated within our own 
halo, but also especially of the component generated during 
the evolution of structure at early times, and visible as an 
isotropic gamma-ray background at the present day. 

In this paper, we calculate the flux produced by an- 
nihilations in simple CDM halos, relative to the flux pro- 
duced in a uniform background. We correct this result for 
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halo substructure, and integrate over a large volume to de- 
termine the cosmological background from WIMP annihi- 
lation as a function of redshift. The outline of this paper 
is as follows. In section 2, we define a dimensionless flux 
multiplier / that accounts for the enhanced rate of two- 
body interactions produced by inhomogeneities in the dark 
matter distribution, and determine its value for simple viri- 
alised halos. In section 3 we calculate / for cosmological 
volumes, using analytic estimates of the halo mass function 
and halo concentrations, and determine its redshift depen- 
dence. Finally, in section 4 we study the contribution to / 
from substructure within virialised halos, and calculate / 
for a set of realistic halos generated using a semi-analytic 
model of halo substructure. Throughout this paper we as- 
sume a Lambda-CDM (LCDM) cosmology with a cosmolog- 
ical constant Ao = 0.7, a matter density fi m ,o = 0.3 and a 
Hubble parameter Ho = h x 100km s~\ with h = 0.65. 



2 ANNIHILATION RATES IN SIMPLE HALOS 

In current hierarchical models, cold dark matter is expected 
to form centrally concentrated halos with a characteristic 
density profile. Dense substructure is abundant within these 
halos, as a relic from earlier stages of the hierarchical merg- 
ing process. Since the annihilation flux is quadratic in the 
density, these inhomogeneities will increase the flux from 
a halo of a given mean density. We begin by computing a 
dimensionless quantity that describes this enhancement. In 
the next section we will then study the evolution of this 
quantity with epoch. 



2.1 The Dimensionless Flux Multiplier 

If dark matter consists of neutralinos with a mass m x and 
a velocity-averaged cross-section for annihilation (erv), then 
the annihilation flux produced within a volume V will be 



<L»xM / f dV . 



(i) 



where p is the local density of CDM. For non-relativistic 
particles, (av) is approximately independent of v so the flux 
will just be proportional to p 2 . 

Since the rate depends quadratically on the density, the 
total rate from a given mass within a given volume will be 
higher if the dark matter is distributed inhomogeneously. We 
can study this enhancement by defining the dimensionless 
flux multiplier 



/GO 



p 2 V 



p 2 dV, 



(2) 



for a distribution within a volume V, where p is the average 
density within this volume. This function can also be written 
as a mass-weighted density average: 



(3) 



where M is the total mass within V. Clearly f — 1 for 
a homogeneous distribution, while for a power-law density 
profile p oc r~ a , 

(3 -a) 2 



f 



3(3 - 2a) 



(4) 



provided a < 1.5. For a = 1.5, integrating equation (y) from 
r m in to r max gives: 



/ = j ln(r max /r mill ) , 

so / diverges logarithmically as r m i n goes to zero. 



(•») 



2.2 Analytic Profiles 

Two analytic density profiles are commonly used to fit the 
spherically averaged properties of dark matter halos, the 
NFW profile (Navarro, Frenk & White 1996, 1997), and 
the Moore profile (Moore et al. 1998). We can specify these 
generically as 

p(r) = p.r? +/h /r a (r fi + rf)"', (6) 

with a = (5 = 1, 7 = 2 for the NFW profile and a = (3 = 
1.5, 7=1 for the Moore profile. The total mass within radius 
r is 



M{< r) = 4nr s p s m(r/r a ) , 
with 

m(x) — ln(l + x) — x/(l + x) 
for the NFW profile and 
m(x) = 2/3 ln(l + a; 1 ' 5 ) 



(7) 
(8) 
(9) 



for the Moore profile, while the mean density within this 
radius is simply 



p(x) = 3p s m(x)/x , 



(10) 



where x = r/r s . 

In a cosmological setting, halos are virialised out to a 
radius r v corresponding to an overdensity A c of roughly 200 
relative to the background. The concentration c = r v /r s 
of a halo describes the size of this radius relative to r a . ^ 
Calculating the flux multiplier over the virialised region of 
a halo, we get a function which depends only on c: 



f(c) 



p 2 V 



p(r) 4nr dr 



3m 2 (c) J x 2a (1 + xP) 2 i 
Thus for the NFW profile 

dx 



x 2 dx. 



f(c) = 



3m 2 (c) 



{1 + xY 
c 3 /(l + c) 3 



(11) 

(12) 

(13) 
(14) 



9(ln(l + c) -c/(l + c)) 2 ' 

with a limiting value of /(c) = 4/3 for small c, whereas for 
the Moore profile, 

.3 

/(c) 



3m 2 ( 



dx 



xil + x 1 - 5 ) 2 



2c" 



9m 2 (c) 



1 + x 1 



(15) 
(16) 



Note that the concentrations derived by fitting the outer re- 
gions of simulated halos differ according to the profile used, with 

C NFW — l-73CMoore- 
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This expression diverges logarithmically as a: m in goes to zero, 
reflecting the fact that the flux from a pure r -1 ' 5 cusp is in- 
finite. In practice, annihilation would reduce the density at 
the centre of the halo to a finite value even for pure CDM, 
while in real halos baryons will also affect halo structure 
on small scales. Furthermore, the simulations that indicate 
central cusps in halos can only probe down to r/r v ~ 0.01, 
or r/r s ~ 0.1 (Power et al. 2002); below this scale, it is 
possible that the slope of the profile continues to change 
(Taylor & Navarro 2001). Thus in what follows we will con- 
sider Moore profiles truncated at various inner cutoff radii, 
either s m i n = 10 -10 (profile Ml hereafter), :r m m = 10 -2 
(profile M2 hereafter), s m i n = 10" 1 (profile M3 hereafter), 
or a; m i n = 1/3 (profile M4 hereafter), with cores of constant 
density within the cutoff radius. (Adding a core modifies 
equations ( |l^ ) and ( |lr^ ) slightly, but the difference is negli- 
gible for x 3> lEmin-) The first cutoff corresponds the scale on 
which annihilation alone might truncate a pure power-law 
density profile (Calcaneo-Roldan & Moore 2000) and the 
third cutoff corresponds to the largest core consistent with 
the results of current high-resolution simulations, while the 
fourth cutoff corresponds to the limit of the region within 
which baryons might plausibly have flattened out the core 
in the halos of massive galaxies. 

2.3 Non-analytic Profiles 

As mentioned above, the simulations that provide evidence 
for a universal density profile for CDM halos can only reli- 
ably resolve two decades of radius below the virial radius. 
While the logarithmic slope appears to converge to a value 
between —1.5 and —1 in the inner regions, it is clearly not 
a simple power-law in the outer regions, where it flattens 
substantially over a decade in radius. Thus, there remains 
the possibility that the inner slope of CDM halos, were it 
resolved to smaller radii, would continue to change in value. 
A shallower inner slope would help to explain the rotation 
curves of dark-matter dominated dwarf galaxies (e.g. Blais- 
Ouellette, Amram, & Carignan 2001; de Blok & Bosnia 
2002), as well as the mass distribution of the Milky Way 
within the solar circle (Binney & Evans 2001). 

One possible form for such a profile was proposed by 
Taylor and Navarro (2001, TN hereafter). In the outer re- 
gions it resembles the NFW or Moore profiles, while in the 
inner regions (below about log(r/r„) = —1.5), its logarith- 
mic slope decreases slowly to —0.75. An analytic fit to the 
inner slope is: 

dlnp _ 0.75 + 2.625 x 1/2 

d \nx ~ 1.0 + 0.5 x 1 ^ ' ^ > 

where x = v/vtn, with rxN = 5/3 r s ^pw- We will calcu- 
late /(c), which is convergent and well-behaved for this pro- 
file, from a direct numerical integration, although it can be 
computed to reasonable accuracy using the fitting formula 
above. From equation (^]), the limiting value for small c is 
/(c) = 9/8. Note that in its outer regions, the equilibrium 
density profile derived by TN goes to zero at a finite radius. 
This assumes an equilibrium that will never be achieved in 
cosmological setting, however, where halos will continue to 
accrete material onto their outer parts. Thus, we assume 
that beyond itn = 1, the profile drops off as r~ 3 , as in the 
NFW and Moore fits. 



2.4 Comparison Between Profiles 

Figure [l] shows the density profiles (top panel) and f(x) 
(bottom panel) for the six profiles discussed above, in- 
cluding four versions of the Moore profile with different 
constant-density cores (within radii indicated by the hori- 
zontal dashed lines). In order to compare the different pro- 
files, they have been plotted in terms of x v = r/r p , the 
radius relative to the point r p at which the circular velocity 
peaks, and normalised so that the density is the same at r p . 

In the top panel, we see that all the density profiles aside 
from M4 are in excellent agreement for x p > 0.1, the re- 
gion probed by current simulations. Existing simulations are 
roughly compatible with any of the fits shown here. The dif- 
ferences in the inferred central behaviour of halos are based 
on different interpretations of the density profile very close 
to the current numerical resolution limit. 

The bottom panel shows the effect of this uncertainty on 
/(c), where c is measured relative to r p , as in the top panel. 
We see that in each case, the flux multiplier increases with 
concentration, reflecting the contribution from the dense 
central core. To first approximation, the dependence on c 
is independent of profile, while the inner slope of the halo 
profile determines the overall normalisation of /(c). For the 
Moore profile it is 2-25 times larger than for an NFW pro- 
file, while for the shallower TN profile it is comparable to a 
Moore profile with a core between 1% and 10% of r s in size. 
The solid vertical lines on the plot indicate typical values of c 
for present-day galaxies and clusters. For galaxy halos, /(c) 
is typically 50 for an NFW profile, 200 for the TN profile 
or profile M2, and 1200 for the most extreme Moore profile, 
Ml. For all profiles it scales with concentration as roughly 
c 2 ' 4 in this range of c, as indicated by the solid line in the 
corner of the figure. 

This figure indicates the importance of the inner region 
in determining the total flux from annihilating dark matter. 
If recent simulations are correct, it seems likely that this 
flux may be 2-4 times larger than that predicted for a pure 
NFW profile; in the extreme case of a Moore profile with an 
inner cusp limited only by self-annihilation, the total flux 
is 25 times the NFW value. Since the dependence on r/r p 
is similar for all profiles, in what follows we will assume 
the NFW form for /(c), and allow that the flux may be a 
constant 2-25 times greater than this. 



3 AVERAGE OVER COSMOLOGICAL 
VOLUMES 

3.1 Local Contribution 

To calculate the background flux from annihilations within a 
large volume, we need to add the relative contributions from 
regions of different density within this volume. Following the 
Press-Schechter approximation (Press & Schechter 1974), we 
can consider all the mass in a given (physical) volume V to 
be contained in virialised halos of some (possibly very small) 
mass. We make the further approximation that these halos 
are spherical, with a fixed virial overdensity A c relative to 
the critical density, and have universal density profiles with a 
concentration c = r v /r s that depends only on their mass. If 
the volume contains V (dn(M)/dM) AAI halos in the mass 
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Figure 1. The dimensionless flux multiplier / as a function of 
concentration c, for the various density profiles shown in the top 
panel. The (solid) vertical lines indicate typical concentrations for 
galaxies and clusters at the present day, while the dashed lines 
indicate the size of the constant-density core in profiles M2, M3 
and M4. See text for a description of the profiles. A line of slope 
2.4 is shown in the bottom right-hand corner for comparison. 



range M to M + AM, then from equation (pj), the total flux 
multiplier for the volume will be 



f(V) 



p 2 V 



V^p- dM ) V/!.. (IS, 



where p is the mean density of dark matter within V, p — 
A c p c is the average density of bound halos, and f(c(M)) is 
the flux multiplier for halos of concentration c = c(M) . Since 
p = fl p c , we can rewrite the dimensionless flux multiplier 
for large volumes as: 



Ac 



(19) 
(20) 



where F(M) is the fraction of the universe in virialised halos 
of mass M or larger. 

F can be estimated using the Press- Schechter formal- 
ism: 



dF P s(M) 
dM 



dM = 



dv 
2x i/2 

7T 



dv. 



(!) 



exp 



dv 



(21) 
(22) 



with v = S c /a(M), where S c is the critical overdensity and 
a(M) describes the power spectrum of density fluctuations. 
Recent simulations (e.g. Jenkins et al. 1998) have suggested, 
however, that this mass function may overestimate the num- 
ber of halos near M», the characteristic mass for which 
a(M) = D(z) 5 C (where D(z) is the linear growth factor 
at redshift z). An alternative mass function, proposed by 
Sheth and Tormen (1999) and based on the ellipsoidal col- 
lapse model, is: 



dFsT{v) 
dv 



dv = 2A 



dFpsiy) 
dv 



dv, 



(23) 



with A = 0.3222 and q = 0.3. 

Since both these mass functions have a power-law be- 
haviour at low masses, where concentrations and thus /(c) 
are larger, it is not clear the integral in equation ( |2o[ ) will 
converge as we include the contribution from smaller and 
smaller halos. (Although the total mass per unit volume 
must converge, / within a given volume need not converge, 



as in the case of a pure 



density profile.) In practice, 



we expect baryonic phenomena to complicate structure for- 
mation on small mass scales, and annihilation itself will also 
limit the contribution from very dense material. Thus, we 
will truncate the integral at some limiting mass, and con- 
sider the behaviour of / as a function of this mass limit. 

The concentration of a halo should reflect the density 
of the universe at the time when it assembled the material 
now in its central core. There are several predictions of the 
concentration-mass relation (Navarro, Frenk & White 1997; 
Bullock et al. 2001; Eke, Navarro, & Steinmetz 2001; Wech- 
sler et al. 2002) , based on this interpretation. We will use the 
most recent analytic model, that of Eke, Navarro, & Stein- 
metz (2001 - ENS hereafter), to calculate the dimensionless 
flux multiplier, though we expect similar results from the 
other models. We use the public code supplied by the au- 
thors to calculate concentrations, and integrate the expres- 
sions above numerically. Note that one point of uncertainty 
in all these models is what minimum concentration to assign 
to halos that have just formed. The ENS code, for instance, 
predicts concentrations of 1 or less for very massive or very 
high-redshift halos, but the least concentrated halos actually 
found in their simulations have c ~ 2-3. For recently formed 
objects the universal profile may not be well established, so 
any uncertainty produced by varying the minimum concen- 
tration reflects the limitations of the analytic model. 

The top panel of figure 2 shows the dimensionless flux 
multiplier contributed by halos over some limiting mass, at 
various epochs in a LCDM cosmology, plotted as a function 
of the mass limit. The solid lines show the results for the 
PS mass function, while the dashed lines show the results 
for the ST mass function. We see that in either case, the 
total flux is more or less convergent as we include smaller 
and smaller masses; at z = 0, it only varies by a factor of 2, 
for instance, for a limiting mass anywhere between 10 10 Mq 
and 10 3 Mq . The value is also reasonably independent of the 
mass function used; it is at most 2 times larger for the PS 
mass function. The uncertainties in halo concentrations are 
only important at high redshift, when many halos have re- 
cently formed. We have actually plotted two sets of curves in 
figure 2, one with a minimum concentration of 1 (the lower 
set of lines at each redshift) and one with a minimum con- 
centration of 2 (the upper set), but the difference between 
them is negligible for z < 10. Finally, the total flux multi- 
plier increases with time, reflecting the progressive growth of 
structure. Of course the density of the universe is decreasing 
at the same time, so the flux relative to that at the present 
day is larger at high redshift; this is shown in the bottom 
panel, where we have used the present-day density to nor- 
malise /. This increased flux will partly compensate for the 
decreased volume element at high redshift, as discussed be- 
low. 
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log[M mjn /M ] 

Figure 2. (Top panel) The dimensionlcss flux multiplier for large 
volumes f(V) as a function of the minimum halo mass considered, 
at different redshifts. The line styles indicate the mass functions 
assumed, while the upper and lower lines visible at high redshift 
indicate the effect of uncertainty in the halo concentrations for 
recently formed objects. (Bottom panel) The flux multiplier rela- 
tive to the present-day density of the universe. NFW halo profiles 
and a LCDM cosmology are assumed. 



3.2 The Universal Flux Multiplier 

In previous section we have derived the flux multiplier for a 
local volume, large enough to contain a representative sam- 
ple of halos, but small enough that evolutionary effects are 
negligible. In a similar fashion we can extend the definition 
of / to entire observable universe. Consider the expression 
for / integrated over volume elements dV(a) at different 
scale factors a, to some maximum a m 



/univ (^m ) 



V(a m ) p(0) 2 



V(a m )p(0) 2 



p(afdV(a) 



p(aff(a)dV(a) 



(24) 



(25) 



where we have chosen to normalise by the present-day den- 
sity of the universe, p(0), /(a) is given by equation (po|), 
and the last term corrects for the time dilation at a given 
epoch. Since a = (l + z)' 1 , this is equivalent to the following 
integral over redshift: 



y*univ(^in) — 

with 



3(l + z m ) 3 r m p(z) 2 f(z)dz 
hi 



(26) 



p(0) 2 J (1 + *) 4 
= (l + 2)- 1 (n m>o (l + 0) 3 +A o )- 1/2 , (27) 



for the LCDM cosmology considered here. 

One quantity of interest is the relative contribution to 
the flux per interval of time, df un i v /dt = (dfumv/dz) \dz/dt\, 
which is plotted as a function of redshift in figure ^. The four 
sets of lines indicate the contribution from mass functions 
truncated at the different minimum masses indicated. The 
line styles indicate the results for PS and ST mass functions. 
We see that the relative contribution peaks at quite early 
redshifts, roughly the epoch when the smallest halos con- 
sidered have formed, as indicated by the arrows. At higher 




Figure 3. The relative contribution to the flux multiplier, per 
interval of time. The four sets of lines indicate the contribution 
from mass functions extending down to the minimum masses in- 
dicated. The line styles indicate the results for PS and ST mass 
functions. A LCDM cosmology is assumed. Arrows indicate the 
redshift at which M* = M m i n , i.e. <r(M m i n ) = D(z)S c . 



redshift, it drops off sharply, because there are very few ha- 
los over the mass limit considered, while at lower redshift 
it drops off more gradually, as small, dense systems merge 
to form larger systems of lower density. The latter drop- 
off may be overestimated, however, since the cores of dense 
halos may not be disrupted by merging, but may survive 
instead as distinct substructure within larger systems. We 
consider this effect in the next section. 



4 ANNIHILATION RATES IN REALISTIC 
HALOS 

In cosmologies dominated by CDM, large halos assemble 
hierarchically through successive mergers between smaller 
objects. This merging process is inefficient, and leaves many 
undigested remnants, or subhalos, orbiting within a galaxy 
or cluster halo. The properties of these objects have been 
studied extensively in several recent, high-resolution SCDM 
and LCDM simulations (Ghigna et al. 1998; Klypin et al. 
1999; Moore et al. 1999; Springel et al. 2001), and can also 
be understood using analytic and semi-analytic models (Bul- 
lock, Kravtsov, & Weinberg 2000; Taylor 2001; Benson et 
al. 2002; Somerville 2002). The latter have the advantage 
that the properties of the underlying halos can be specified 
directly, and are independent of resolution and numerical 
effects. In this section we use a recently developed semi- 
analytic model of halo formation (Taylor & Babul in prepa- 
ration; see also Taylor 2001) to determine the flux multiplier 
for a realistic halo with substructure, at various different 
epochs. 



4.1 The Contribution from Halo Substructure 

From equation ^, the relative contribution to the flux multi- 
plier from an individual subhalo within a larger halo is pro- 
portional to p s hM s h/'(cah)/pbgMbg/(cbg), where the sub- 
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script sh indicates the subhalo and the subscript bg indi- 
cates the average for the background system in which it re- 
sides. This implies that if subhalos are substantially denser 
or more concentrated than the main system, they can domi- 
nate /, even though they only contribute to a small fraction 
of the total mass of the system. This is particularly true 
since tidal stripping will remove the lowest-density material 
from substructure, increasing the relative contribution from 
the material which remains bound. 

More precisely, if we decompose the density of CDM 
within a halo into two components, a smooth background 
following the basic profiles given in section 2, and a set of 
subhalos which have been stripped to varying degrees, then 

p(x) = Pbg(x) + ^ Psh,i( x )- ( 28 ) 

i 

The contribution to / then consists of several terms 

/ = -Ar i fo* + Y, p^ 2 dV ( 29 ) 

/ (pbg + Pb S Psh,i + dv (30) 



]>>w Xi )+]T^-/'( Cl 



(31) 



where Mb g , Cb g and Vb g are the mass, concentration and 
volume of the main system, and Mi, a and Vi are the cor- 
responding properties of each subhalo. In deriving the final 
result, we have assumed that the subhalos do not overlap, 
and that each one is small compared to the main system, so 
that pbg is more or less constant within the volume covered 
by a given subhalo. The function /'(c) is the dimensionless 
flux multiplier defined in section 1, modified to account for 
the effects of tidal heating and stripping, which may change 
the density profile of the subhalos as they orbit within the 
main system. We derive an expression for this function in 
terms of the original profile and the amount of mass lost 
below. 



4.2 The Flux Multiplier for Subhalos 

The density profiles considered in section (Q) characterise 
isolated halos; interactions between halos or mergers with 
larger systems will modify these profiles through a combi- 
nation of tidal heating and tidal mass loss. Numerical sim- 
ulations show that tidal effects tend to strip mass off the 
outer regions of a halo preferentially, preserving the central 
density of the halo until it has lost a large fraction of its 
original mass. Hayashi et al. (2002) propose the following 
modification of the density profile to describe this process: 



p'(r) 



ft 



1 + (r/r t ) 3 



p(r) 



(32) 



where rt is the tidal radius to which the halo is stripped, and 
ft is the reduction in the central density. For a halo with an 
NFW profile and an original concentration of 10, they find 
that these parameters are related to the fraction rrih of the 
original mass still bound to the halo by 



log (r t ) 
and 



1.02 + 1.38 log(ro b ) + 0.37(log(m b )) 



(33) 



1 1 1 

V = V(< 10 r a ) 




- V = V(< r L ) 




V = V(< 3 r L ) 

1 





Figure 4. The dimensionless flux multiplier /' calculated over 
various volumes, for stripped systems with an original NFW pro- 
file of concentration cq = 10, as a function of their stripped con- 
centration ct = rt/r 3 . The dashed lines show /' calculated within 
the original tidal radius rt,o = 10 r a , the dotted line shows /' for 
the volume within rt, and the dot-dashed line shows /' for the 
volume within 2 rt . The solid line shows / for the original density 
profile, truncated at different radii. 



log(/ t ) ~ -0.007 + 0.35 log(m b ) + 0.39(log(m b )) 5 
+0.23(log(m b )) 3 . 



(34) 



Using this modified density profile, we can calculate /', 
the flux multiplier for tidally stripped halos. This should 
be a function of Co, the original concentration of the halo, 
and Ct = r t /r s , the concentration of the stripped system. 
In practice, we shall consider /'(ct;co = 10) since this was 
the profile modeled by Hayashi et al. (2002). In defining /', 
we also need to choose an appropriate volume over which 
to integrate. Since the modified density profile in equation 
( [32] ) is not truncated completely at rt, we could continue to 
integrate over the original volume; in this case, /'(c) would 
rise sharply (roughly as (rt/r^o) 3 ) as the mass of the halo 
is confined to a smaller and smaller region of the original 
volume. Since rt roughly defines the region in which the 
subhalo density profile dominates over the background, the 
approximation made in deriving equation ([u]) will be more 
accurate if we define /' over the volume within rt- In any 
case, in this section we will show the results of the calcu- 
lation for several possible choices of volume, to clarify the 
effect on /'. 

Figure 4 shows f'(c t ;co = 10) for the volume within 
r/r s = Co = 10 r s (dashed line), r — rt (dotted line), and 
r = 2r t (dot-dashed line), with f(ct) shown for comparison 
(solid line). The truncated functions differ from the origi- 
nal one at Ct = Co = 10 because the density profile given by 
equation ([32]) is different from the unstripped one even when 
r t /r s — cq. We see however that if we consider the volume 
within r t (dotted curve), /'(c t ;co = 10) ~ 1.3 f(ct) to good 
approximation. We will assume that this result remains true 
for halos with larger or smaller original concentrations, pro- 
vided Ct < 0.9 Co, and use this to calculate the contribution 
to the total flux from evolved subhalos. 
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Figure 5. Top panel: the contribution to the flux multiplier con- 
tributed by substructure over a minimum mass M m j n , as a func- 
tion of M m j n , relative to the flux multiplier f\, E expected for a 
smooth NFW halo without substructure. The solid curves are av- 
erage results for a set of semi-analytic galaxy halos with a total 
mass of 1.6 xlO 12 Mq, in a LCDM cosmology. The dashed vertical 
line indicates the lower mass limit of the merger trees used. Bot- 
tom panel : as above, but showing the redshift dependence of the 
relative contribution for two specific mass limits, Mjj m = 10 7 Mq 
and M lim = 10 s Mq. 



4.3 Results for LCDM Halos 

Using equations ( pi] ) - (pS4|), we can calculate the total flux 
multiplier for composite systems with substructure. To do 
so, we use a semi-analytic model of halo formation (Taylor & 
Babul in preparation; see also Taylor 2001) to determine the 
number, masses, density profiles and positions of subhalos 
within a set of hierarchically formed halos. We have followed 
halo substructure down to ~ 10~ 5 times the mass of the 
main system. Since the spectrum of substructure is fairly 
steep at low masses, the total flux multiplier may depend 
on this cutoff, so we will consider our results as a function 
mass limit, as in section (|^), and extrapolate downwards to 
estimate the total flux multiplier for smaller limiting masses. 

The top panel of figure |B| shows the relative contribu- 
tion to the flux multiplier contributed by substructure over a 
minimum mass limit M m i n , as a function of M m i n , in models 
using NFW profiles. The contribution has been normalised 
to the value /b g expected for a smooth NFW halo with- 
out substructure. The solid curves are average results for 
a set of semi-analytic galaxy halos with a total mass of 
1.6 x 10 12 M Q , in a LCDM cosmology. The dashed vertical 
line indicates the resolution limit of the merger trees used; 
below this, the spectrum of subhalos becomes increasingly 
incomplete, particularly below 2 x 10 7 Mq. At high redshift, 
where the halos are forming through major mergers around 
this mass scale, the curves are steep so the flux estimate will 
depend sensitively on the mass limit. At low redshift, the ef- 
fect of the mass limit is less important, as most of the flux 
comes from more massive subhalos. In either case, however, 
the total flux appears to have a power-law dependence on 
Mii m at low masses, so it is straightforward to extrapolate 
our estimates to lower mass limits. 

Overall, we see that substructure makes an important 




Figure 6. As figure but for two specific mass cutoffs (solid 
lines), and including the effects of substructure on the total flux 
multiplier (dashed lines). The arrows indicate the redshift at 
which M t = M m i n , as before. 



contribution to total flux right up to the present day, when it 
exceeds the contribution from the background halo. The bot- 
tom panel shows the redshift dependence of the contribution 
more clearly, for two different mass limits, M\i m = 10 7 Mq 
and M\i m — 10 8 Mq . As expected, at high redshift the con- 
tribution from substructure increases somewhat faster as we 
lower Mn m , reflecting the steep slope of the high-redshift 
curves in the top panel, while at low redshift the increase is 
slower. If we extrapolate the power-law dependence seen in 
the top panel to smaller values of M\\ m such as Mn m ~ 10 4 - 
10 , we expect a net contribution from substructure of any- 
where from 10 to 20 times the smooth background value for 
the NFW profile. 

We can determine the effect of substructure on the 
gamma-ray background by cor recting the relative contribu- 
tion calculated in section (3.2) for the additional contribu- 
tion from substructure in shown in figure || this is indicated 
by the dashed lines in figure 0, for halos with NFW profiles, 
and mass cutoffs of Mu m = 1CT Mq and Mu m = 10 s Mq. We 
see that when substructure is taken into account, the flux 
multiplier still peaks just after the epoch when M* = Mii m 
and declines at lower redshift, but both the peak value and 
the value at z = are substantially enhanced. We note that 
the peak and subsequent decline may partly compensate 
for the cosmological factor \dt/dz\ and for intervening ab- 
sorption, which reduce the high-redshift contribution to the 
gamma-ray background. We will investigate the detailed line 
and continuum spectra predicted by this model in a subse- 
quent paper. 



5 DISCUSSION 

We have calculated the amount by which structure in dark 
matter on subhalo, halo and cosmological scales amplifies 
the expected signal from neutralino annihilation. The overall 
magnitude of the effect remains uncertain at several levels, 
some more important than others. Uncertainty in the fun- 
damental density profile which characterises halos, or more 
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specifically in the slope of their density profiles at small radii, 
can change the flux multiplier for simple halos by a factor of 
20-25, assuming concentrations typical of galaxy halos (see 
figure |l|) . There is evidence in recent simulations for some 
excess mass in the inner regions, relative to the NFW pro- 
file (Power et al. 2002), so a reasonable estimate for the flux 
multiplier may be 3-4 times the NFW value used for most of 
the calculations in this paper, although studies of the Milky 
Way (e.g. Binney & Evans 2001) suggest less dark matter 
in its central regions than these profiles would imply. There 
is no direct evidence for a pure Moore profile extending to 
very small scales, although if this is the case the resulting 
flux from simple halos will be 20-25 times the NFW value 
used here. 

The mass function of dark matter halos, averaged over 
large volumes, is also uncertain, particularly at the low mass 
end. The mass function determined from simulations is bet- 
ter fit by the form proposed by Sheth and Tormen (ST) 
than by the traditional Press and Schechter (PS) mass func- 
tion. Using the ST mass function to predict the background 
flux reduces its amplitude by about 1.5 over PS. We assume 
that the scale invariance of dark matter substructure will be 
broken at very small masses, perhaps close to Jeans mass 
at recombination, for instance. At a minimum, there should 
be some maximum density for dark matter halos, or equiv- 
alently some maximum for a(M) at early times. Varying 
the lower mass limits over a reasonable range (10 4 -10 8 Mq) 
changes the flux by a factor of 2 or so for simple halos, but 
the contribution from substructure may increase this differ- 
ence to as much as a factor of 40. 

Overall, these uncertainties combine to produce a to- 
tal uncertainty of more than an order of magnitude in the 
flux multiplier. A conservative model, at the bottom of this 
range, consists of an NFW profile, with a ST mass function 
and a mass limit of 10 6 Mq . For this model, the present-day 
flux multiplier is 6 x 10 5 . We estimate that a more likely 
model has a profile with more mass in its inner regions, and 
lower limit to the mass function of 10 5 Mq. For this model, 
the present-day flux multiplier is 5 x 10 6 . For an extreme 
model, with a Moore profile limited only by annihilation, 
and a mass limit close to the Jeans mass at recombination, 
the present-day flux multiplier will be ~ 1 x 10 8 . 

It is worth noting that with these large flux multipliers, 
some neutralino candidates can already be ruled out. The 
86 GeV neutralino considered in Bergstrom, Edsjo, & Ullio 
(2001), for instance, produced almost 10% of the gamma-ray 
background observed by EGRET at 1 GeV, assuming a flux 
multiplier of ~ 2 x 10 6 , while a 166 GeV neutralino produced 
more than 50% of the flux observed at 10 GeV. Both these 
candidates are thus possible in our most conservative model, 
while the higher-energy candidate is marginally excluded in 
our favoured model, and both are ruled out in our extreme 
model. 
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